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Abstract 

Recently there has been a claim on the complete irrelevance of quan- 
tum modeling for plasmas. We address this subject from basic prin- 
ciples. Physical situations where quantum effects play a decisive role 
are identified. 
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Recently Vranjes, Pandey and Poedts wrote a series of comments [U [2] 
on the limited applicability of quantum effects in plasmas, except for atomic 
and/or nuclear aspects like ionization, calculation of scattering cross sections 
and so on. In spite of the prompt reply by Shukla and Akbari-Moghanjoughi 
[3], besides the existing thorough discussions on the basic theoretical models 
for quantum phenomena in plasmas [H El El [7] , we fell it would be still of some 
interest to have more independent views on the topic. Hence we address the 
question from first principles, apologizing for the superposition with some 
well-known arguments. 

For any collection of particles, one may roughly separate quantum effects 
in three categories: (a) the class of quantum effects due to particle indis- 
tinguishability; (b) the family of phenomena originating from wave-particle 
duality; (c) the class of quantum effects linked to the intrinsic magnetic 



moment of the charge carrier dynamics, irrespective of quantum statistical 
issues. Let us start discussing point (a), with an eye to systems of massive 
fermions since usually electrons are an essential component of plasmas. In 
this case, there are two main possible descriptions which are in order, based 
resp. on classical, Maxwell-Boltzmann (MB) or quantum, Fermi-Dirac (FD) 
statistics. The MB is an approximation to the FD statistics, valid when there 
is not too much wave packet overlapping. The condition for the applicability 
of MB is given by T 3> Tp, where T is the thermodynamic temperature and 

T F = -^(3vr 2 n ) 2 / 3 (1) 
2mK B 

is the Fermi temperature. In Eq. (jTJ), K is Planck's constant divided by 
2tt, m is the charge carriers mass (in this text we consider electrons), Kb 
is Boltzmann's constant and n is a reference equilibrium number density. 
Apart from an irrelevant numerical factor of order unity, from Eq. ([I]) the 
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sufficiently high temperature condition is found to be equivalent ton ^> 
Xb, where 

B \f2TxmKB~T ^ ^ 

is the thermal de Broglie length for a Fermi gas. Therefore, if the mean 
inter-particle distance is much larger than the typical wave packet size one 
may disregard the indistinguishability aspects. For simplicity, we are talking 
about a completely ionized gas; for bound particles some appropriate length 
scale should replace the thermal de Broglie length as a measure of the wave 
functions size. 

In conclusion, a many-body fermion system deserves a quantum statistics 
treatment whenever it is sufficiently cold or dense. Accordingly we refer to de- 
generate (FD) and non-degenerate (MB) plasmas. We don't intend to present 
an exhaustive list of real degenerate plasmas, limiting ourselves to few exam- 
ples: (i) the plasma appearing in intense laser compression schemes jH], which 
nowadays can have large densities of the order no — 10 32 m -3 . In this case one 
need the unlikely elevated temperature T ^> 10 7 K to be allowed to apply MB 
statistic without committing serious error. One should note that with the in- 
creasing laser power entering the multi-petawatt range one can even envisage 
to probe strong quantum aspects of nature, including nonlinear modifications 
of the Maxwell equations due to photon-photon scattering mechanisms. Such 
avenues, however, can be accessed only through electric field strengths of the 
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order of the Schwinger limit m 2 c 3 /(eft) ~ 10 18 V/m, where c is the speed of 
light and e the elementary charge. Quantum field theoretical aspects are far 
outside the scope of most quantum plasma models; (ii) the plasma in com- 
pact astrophysical objects like white dwarfs [H] and neutron stars [TU] where 
n ~ 10 36 m~ 3 . In this case FD statistics is necessary unless T 10 9 K, 
which is also unlikely. 

At least one can briefly mention some further quantum plasma examples, 
one of the most notorious being the degenerate electron gas in metals as re- 
membered in [3J. As remarked in [2J, the properties of such quantum plasmas 
have been already discussed decades ago. Also one can consider degenerate 
plasmas in the interior of giant planets like Jupiter, plasmas confined in traps 
at sub-kelvin temperatures, and ultra-cold plasmas (with T < lO -3 ^) gener- 
ated from Rydberg states. We urge the reader to examination of the Figure 1 
of Refs. [5J, [3 [TTJ d2] , showing density versus temperature diagrams where the 
MB and FD regions and specific physical systems in both classes are clearly 
identified. Notice that in pQ a similar density versus temperature is shown, 
but restricted to densities of the order of 10 17 m~ 3 , which is obviously not 
representative as already pointed out in [3J. Following the comments in [7J, 
from more complete diagrams it is easy to realize that the family of Coulomb 
systems (many-body systems where the electric force plays a dominating 
role) has grown far beyond the conventional laboratory and space plasmas. 
In the same sense, the very concept of "plasma" is continuously evolving. 
In this context, the relevance (or irrelevance) of quantum statistical effects 
needs by-case analysis in each concrete situation. 

Let us now turn our attention to the second class of quantum effects 
which can impart the behavior of a many-body system, namely class (b) 
due to wave-particle duality. This includes the whole manifold of quantum 
diffraction effects, appearing already in the framework of the Schrodinger 
equation (in the non-relativistic case), without one-to-one relation with the 
Pauli exclusion principle. For instance, wave function spreading and tun- 
neling are among the eminent undulatory processes which are of recognized 
relevance in the dynamics of ultra-small semiconductor devices. In this res- 
pect, consider the Wigner-Poisson model [T3] for electrostatic fields in one 
spatial dimension, 

d£ p_d£ ie^ fdv'dcazv( i(p ~ P ' )s ) x 
dt + mdx + 2nh 2 J dpdSG V { h ) X 
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x (x+^-(f> (ac-f /(^,P'^) = 0, (3) 

= ^Jdp'f(x,p',t)-n (xfj . (4) 

In Eqs. flS])-®, / = f(x,p,t) is the Wigner function, = <f)(x,t) is the 
scalar potential, e is the permittivity and n (x) is a position-density ionic 
background, to be specified according to the doping profile. In spite of at- 
taining negative values in certain regions in phase space (since Eq. ([3]) does 
not preserve the positive definiteness of /), the Wigner function can be used 
to evaluate macroscopic quantities like number, current and energy densities 
basically in the same manner as a classical probability distribution function. 
It is well-known [T3j that the Wigner-Poisson system is a convenient tool to 
describe the resonant tunneling diode (RTD). As the name says, the RTD 
relies on the quantum effect of tunneling, hence no classical model like the 
Vlasov-Poisson one could be appropriate for it. This is independent of the 
carrier concentration, although the Fermi-Dirac character needs to be taken 
into account in the drain region of n + nn + diodes for instance. In applica- 
tions, we have both FD and MB statistics when dealing with such nanoscopic 
systems [13J. 

Whether one call a charged particle system described by Eqs. ©-AID a 
"plasma" or "an electron gas in a semiconductor" is just a matter of taste. 
At the end the basic equations are the same, perhaps with the difference that 
usually macroscopic plasmas consider an homogeneous ionic background. Ac- 
tually the position-dependence of the doping rio(x) in Eq. (T4]) makes the anal- 
ysis nonlinear ab initio, in contrast to the homogeneous case where one can 
start searching for linear waves and instabilities, just as in classical plasma 
theory. There is no profound reason to discard a complete analysis of the 
Wigner-Poisson and similar quantum kinetic models from the plasma physics 
point of view. Such a deleterious attitude would sound like an artificial cen- 
sorship. Nevertheless, in certain situations [TS] a clever order of magnitude 
analysis justify the neglect of the extra dispersion associated to the nonlocal 
term in Eq. ([3]), while keeping other quantum terms originating e.g. from the 
intrinsic spin dynamics. However, presently we are far from having a com- 
plete picture of quantum plasma phenomena. Therefore a systematic neglect 
of quantum diffraction effects from the very beginning is not advisable. 

Equations ©-dlj) constitute an integro-differential system, to be com- 
plemented with suitable initial and boundary conditions. Moreover, even if 
the Wigner-Poisson system is one of the simplest possible quantum kinetic 
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models (non-relativistic, electrostatic approximation, no spin effects etc.) 
by inspection of Eqs. (J3J)- (J3D it is clearly out of reach to access nonlinear 
phenomena in this framework, except perhaps from numerical simulation. 
This justifies the search for alternative models involving e.g. effective quan- 
tum potentials, time-dependent density functional theory or hydrodynamic 
equations. Each strategy has its merits and limitations. Focusing on hydro- 
dynamic models, they can be deduced taking the moments of the quantum 
Vlasov equation Q, in the same way as the classical fluid equations can be 
derived from the moments of the Vlasov equation [TBj. The fact that the 
quantum Vlasov equation, as it stands, does not have collision terms is im- 
material. The final result would be the same after adding a suitable quantum 
Boltzmann-like collision integral in the right-hand side of Eq. ([3D , as long as 
this collision term is mass, momentum and energy preserving. In this respect 
notice that it is more direct to add phenomenological dissipation terms in 
the fluid equations, if necessary. 

In the same way as in the transition from classical kinetic to classical hy- 
drodynamic models, one faces a closure problem, since the evolution equation 
for the moment of order N involves the moment of order N + 1 of the Wigner 
function. One way to circumvent this difficulty is postulating an equation 
of state, compatible with the linear dispersion relation obtained from kinetic 
theory. However, it is evident that using fluid equations always imply a 
loss of information; the same holds in the classical scenario. Nevertheless, 
perhaps no one could seriously propose the understanding of the nonlinear 
structures in quantum plasmas based solely on numerical simulation of the 
kinetic equations. An interplay of kinetic theory and the analytic insights 
coming from the fluid approach is advisable. Note that quantum hydrody- 
namic models are ubiquitous e.g. in the analysis of ultra-small electronic 
devices and quantum chemistry [T3], seemingly always in the electrostatic 
approximation. 

To make these points more clear, we will explicitly write the quantum 
Vlasov equation for the Wigner function in the presence of electromagnetic 
fields assuming the Coulomb gauge, 

x [0(r+ * *)-0( r /(r,p',t) 
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^^V-JJdsdpe . x 
[A(r+|,f) + A(r-|,f)]/(r,p',f) 



x 



«e f f i i i '(p-p')- s 

- — - — — - p • / / as dp e a x 
hm(2-nhf v J J v 

[A(r+~ t)-A(r-|,t)]/(r, P ',t) = 



(5) 



where f — f(r,p,t),<p — <p(r, p, £) and introducing the vector potential A = 
A(r, t). Equation fl5J) is coupled with Maxwell equations with self-consistent 
charge and current densities determined from the moments of the Wigner 
function. Since magnetic fields are an essential ingredient in plasmas, a 
quantum kinetic theory for plasmas should necessarily include Eq. (EJ) or 
some variant of it. Obviously it is out of reach to progress from this pers- 
pective, except for very simple cases (homogeneous magnetic fields, linear 
theory). Even the numerical simulation of the resulting quantum Vlasov- 
Maxwell system is clearly a challenge. To our knowledge there is no available 
3D code for such a problem, which is still somehow limited (no relativistic nor 
spin effects, no exchange-correlation interactions taken into account). Facing 
this situation, one can: (i) postulate that magnetic fields does not exist in 
quantum plasmas and restrict forever to a peaceful electrostatic world; (ii) 
suggest clever approximations. 

Following the second alternative, one can for example proceed to quan- 
tum hydrodynamic models including magnetic fields. One may question the 
correct form or the applicability domain of such macroscopic formulations 
(as was detailed e.g. in [SI E]); it is however absurd to question the need of 
simplified treatments, after open minded examination of Eq. §5§ or any of 
its equivalents. In an optimistic view, one still could imagine that magnetic 
fields in physically interesting situations could be treated in terms of purely 
classical models, letting quantum diffraction to the electric field alone. How- 
ever, this would imply significant error, as in the case of radio wave dispersion 
in a pulsar magnetosphere [17J, where we have the parameter B/B c of order 
unity, where B is the magnetic field strength and B c = m 2 c 2 / (eh) ~ 10 9 T is 
the critical magnetic field. 
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In close connection to magnetic field effects, one further quantum aspect 
deserving attention is the point (c) raised at the starting of this work: the 
dynamic of intrinsic magnetic moment of the charge carriers, or spin, see [Uj 
for a review. The relevance of the spin dynamics should be evaluated in each 
particular instance. However, in this regard a representative dimensionless 
parameter is given by fisB /(kbT), which is the ratio between the Zeeman 
energy associated to the ambient magnetic field and the thermal energy, 
where \xb = eh/ (2m) is the Bohr magneton. For a sufficiently high magnetic 
field (such as B ~ 10 10 T in the vicinity of a magnetar), or for sufficiently low 
temperature, one can have a significant role of the spin dynamics for instance 
in the propagation of Alfven waves 

To conclude, in this work we have shown that quantum plasma systems 
are ubiquitous in nature and briefly discussed their modeling. It is apparent 
that quantum plasma physics is an exciting, under construction area, in the 
need for creative approaches, both theoretical and experimental. We have not 
discussed the challenge of an efficient treatment of non-ideality (collisional) 
and relativity terms. 
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